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the  results  of  this  study  stimulate  more  research  towards  predicting 


and  improving  the  lift-to-drag  ratio  of  circulation  controlled,  high- 
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Abstract 


Wind  tunnel  tests  were  conducted  to  determine  the  effects  of 


splitter  plate  position  on  trailing  edge  modifications  of  a circulation 


controlled  airfoil.  Analytic  studies  were  also  conducted  to  determine 


the  feasibility  of  using  a potential  flow  computer  program  to  predict 


the  results  of  the  wind  tunnel  tests.  The  airfoil  model  was  elliptical 


in  shape,  20  percent  thick,  and  had  five  percent  camber.  It  employed 


a blowing  slot  for  circulation  control  and  a splitter  plate  for  reduc- 


tion in  mixing  losses.  Modifications  included  slot  positions  on  the 


upper  surface  of  96  and  97  percent  chord,  slot  angles  of  +5  and  -33 


degrees,  circular  and  elliptical  aft  contours,  and  splitter  plate  posi- 


tions on  the  lower  surface  of  99  and  95.  3 percent  chord.  Tests  were 


conducted  at  a Reynolds  number  of  7.4  x 10  and  blowing  momentum. 


coefficients  of  zero  and  0.  03. 


The  lift-to-drag  ratio  of  the  airfoil  was  increased  with  the  use 


of  a splitter  plate  over  a clean  configuration  and  was  maximized  with 


the  splitter  plate  at  the  99  percent  chord  position.  These  .results 


held  true  for  the  variations  in  slot  position,  slot  angle,  aft  contour, 


and  blowing  rate. 


The  potential  flow  computer  program  was  capable  of  matching 


lift  coefficients  obtained  in  wind  tunnel  tests.  However,  due  to  an 


inviscid  flow  assumption,  the  program  was  not  feasible  for  prediction 


purposes  because  it  did  not  model  the  flow  accurately. 


‘ . l . It 'J. V.  _ mA 
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AN  ANALYTIC  AND  EXPERIMENTAL  STUDY  OF  THE 


EFFECTS  OF  SPLITTER  PLATE  POSITION  ON 


THE  TRAILING  EDGE  MODIFICATIONS  OF  A 


CAMBERED  CIRCULATION  CONTROLLED 


ELLIPTICAL  AIRFOIL 


I.  Introduction 


During  the  past  several  years  a great  interest  has  developed  in 
aeronautical  structures  that  improve  low-speed,  high-lift  capabilities. 
The  motivation  for  this  interest  is  the  desire  that  aircraft  land  and 


takeoff  in  shorter  distances  than  normal.  Short-takeoff-and-landing 
aircraft  (STOL)  have  this  as  a requirement  for  their  missions,  as 
they  must  enter  and  leave  areas  normally  inaccessible  to  aircraft. 
Reconnaissance  and  forward  control  aircraft  have  need  of  the  ability 
to  loiter  around  target  areas,  therefore  requiring  a low-speed,  high- 
lift  capability.  Supersonic  fighters,  primarily  designed  for  high 
speed  flight,  normally  do  not  perform  well  at  the  lower  speeds  and 
need  more  lift  during  landing  and  takeoff.  Thus,  high-lift  structures, 
or  wings,  that  perform  well  at  lower  speeds  would  greatly  enhance 
these  missions. 


Poor  performance  at  lower  speeds,  where  higher  lift  is  needed, 
is  due  to  separation  of  the  flow  around  the  wing.  This  is  a result  of 
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the  fact  that  the  flow  loses  energy  near  the  trailing  edge  of  the  wing 
at  high  angles  of  attack  and  can  no  longer  remain  attached  to  the 
surface.  This  separation  of  flow  actually  decreases  the  lift  of  the 
wing,  and  thus  of  the  whole  aircraft. 

Previous  Work 

Through  studies  done  in  the  past,  it  is  known  that  the  flow 
around  a body  can  be  re-energized  to  help  it  remain  attached  to  the 
surface  at  high  angles  of  attack.  The  ways  of  re-energising  this  flow 
before  it  separates  are  called  circulation  control  methods,  and 
include  suction,  blowing,  and  cooling  of  the  surface.  The  blowing 
method  is  the  primary  concern  of  this  present  study.  Blowing,  or 
injecting  air  into  the  already  existing  flow,  adds  energy  to  the  flow'  by 
increasing  its  velocity  and  helps  prevent  it  from  becoming  detached 
from  the  surface. 

This  method  of  circulation  control  was  emphasized  in  previous 
experimental  studies  by  many  researchers.  Kind  and  Maull  investi- 
gated a low-speed,  circulation-controlled,  elliptic  airfoil  with  blow- 
ing near  the  trailing  edge,  and  achieved  fairly  high  coefficients  of 
lift  (Ref  2).  Englar  (Ref  1),  with  another  elliptic  airfoil,  achieved 
higher  lift  coefficients  than  Kind  and  Maull.  Williams  and  Howe, 
with  a slightly  thicker  and  more  cambered  airfoil  than  Englar, 
achieved  even  higher  coefficients  of  lift  (Ref  6).  Stevenson  (Ref  5) 
and  Rhynard  (Ref  4),  at  the  Air  Force  Institute  of  Technology  (AFIT), 

2 


'I*  «6awk.* 


E 


GAE/AE/75D-  12 


made  similar  studies  but  they  also  made  use  of  a splitter  plate.  The 
splitter  plate  was  mounted  on  the  lower  surface  near  the  trailing  edge, 
and  it  spanned  the  model.  It  was  designed  to  reduce  mixing  losses 
near  the  trailing  edge  and  resulted  in  higher  lift  coefficients  and  lower 
drag  coefficients.  Thus,  there  is  much  known  about  circulation  con- 
trol, but  still  remains  a large  area  to  investigate,  even  with  the 
blowing  method. 


Present  Studv 


The  purpose  of  this  study  was  to  further  the  investigation  of 
the  circulation  control  method  of  blowing  air  onto  the  trailing  edge. 
The  principal  goal  was  to  study  the  effec.  of  splitter  plate  positions, 
or  lower  surface  configurations,  on  various  methods  of  flow  accel- 
eration and  trailing  edge  flow  treatment.  An  analytic  study  was  also 
conducted  through  the  use  of  potential  flow  theory  and  computational 
methods  in  an  attempt  to  predict  the  re  alts  of  these  variations. 


Scope 

This  study  was  conducted  in  two  phases.  The  first  phase  con- 
sisted of  researching  and  formulating  feasible  model  modifications 
that  would  produce  variations  in  flow  acceleration  and  trailing  edge 
flow  treatment.  Variations  in  flow  acceleration  were  made  by  chang- 
ing the  blowing  slot  location  and  angle.  Flow  treatment  variations 
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were  made  by  changing  the  aft  contour.  The  splitter  plate  positions 


to  be  used  in  testing  were  then  determined.  This  phase  also  included 


the  development  of  the  potential  flow  model  that  would  be  used  to 


predict  results  of  the  wind  tunnel  testing. 


The  second  phase  of  the  study  was  the  testing  of  the  model  with 


modifications  in  the  AFIT  Five-Foot  Wind  Tunnel.  It  included  the 


manipulation  of  the  potential  flow  computer  program  to  predict  the 


results  achieved  in  the  actual  testing.  The  wind  tunnel  testing  was 


conducted  with  modified  endpieces  that  had  only  one  variation  each, 


from  the  original,  of  the  slot  position,  the  slot  angle,  and  aft  con- 


tour. The  configuration  on  the  trailing  edge  lower  surface  was  also 


varied  for  each  endpiece  from  clean  to  a splitter  plate  in  the  99  per- 


cent chord  position,  and  to  a splitter  plate  in  the  95.  3 percent  chord 


position.  The  tests  were  conducted  through  a range  of  angles  of 


attack  from  -6  to  +6  degrees,  at  blowing  momentum  coefficients  of 


zero  and  0.  03,  and  at  a Reynolds  number  of  7.  4 x 10.  The  vari- 


ations in  the  endpieces,  the  lower  surface,  the  angles  of  attack,  and 


momentum  coefficients  were  incorporated  into  the  potential  flow 


computer  program.  The  resulting  coefficients  of  lift  were  then  com- 


pared to  the  actual  wind  tunnel  data  to  determine  if  the  method  was 


suitable  for  prediction  purposes. 


GAE/AE/75D-12 


II.  Preliminary  Investigation 


The  purpose  of  the  preliminary  investigation  was  to  determine 


what  modifications  were  to  be  made  on  the  airfoil  and  what  lower 


surface  configurations  were  to  be  used  in  the  wind  tunnel  testing. 


Previous  studies  were  investigated  and  it  was  determined  that  vari- 


ations of  three  characteristics  of  the  trailing  edge  of  the  model  should 


be  considered.  These  characteristics  were  the  blowing  slot  position, 


its  angle,  and  the  aft  contour. 


The  slot  location  is  an  important  factor  in  improving  lift- to- 


drag  ratios.  The  further  forward  the  slot  is  located,  the  larger  the 


momentum  coefficient  needed  to  effectively  turn  the  flow  around  the 


trailing  edge.  The  further  aft  the  slot  is  located,  the  larger  the 


momentum  coefficient  needed  to  prevent  separation  of  the  flow  at  high 


angles  of  attack.  Thus,  a variation  in  the  slot  position  was  consid- 


ered and  two  positions  were  chosen. 


The  slot  angle  is  also  an  important  factor  in  improving  1 if t- 


to-drag  ratios.  The  more  tangential  the  air  is  injected  into  the 


surface  flow,  the  more  effectively  it  re-energizes  that  flow  and  helps 


it  remain  attached  to  the  surface  around  the  trailing  edge.  There- 


fore, a change  in  the  slot  angle  was  also  considered  and  two  angles 


were  used. 
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The  effect  of  the  aft  contour  had  been  investigated  previously 
in  the  studies  of  Kind  (Ref  2),  Williams  (Ref  6),  and  Englar  iRef  1). 
They  lound  that  the  turning  of  the  flow  was  more  effective  with  a more 
circular  trailing  edge.  However,  this  more  effective  turning  also 
produced  greater  losses  when  the  upper  and  lower  surface  flows 
mixed.  These  results  led  to  more  lift  out  increased  drag.  Thus,  an 
optimum  combination  of  .'.ircular  and  elliptic  trailing  edge  contours 
was  recommended.  Based  on  this,  two  contours  were  chosen. 


The  last  area  studied  was  the  lower  surface  configuration.  In 
the  investigations  of  Rhynard  (Ref  4)  and  Stevenson  (Ref  5),  a splitter 
plate  was  used  on  the  lower  surface.  The  splitter  plate  was  designed 
co  reduce  the  losses  due  to  mixing  of  the  lower  and  upper  surface  air 
flows  as  they  neared  the  trailing  edge.  Stevenson  determined  the 
necessity  of  such  a splitter  plate  and  Rhynard  attempted  to  optimize 
the  position  on  the  lower  surface  and  the  angle  of  the  splitter  plate. 
Based  on  Rhynard's  results,  three  configurations  were  used  in  this 
study. 

As  a result  of  this  preliminary  investigation,  four  endpieces, 
to  be  mounted  on  the  airfoil,  vere  designed  and  constructed  for  tests 
in  che  AFIT  Five-Foot  Wind  Tunnel.  The  endpieces  differed  from 
each  other  by  a single  variation  in  one  of  the  characteristics.  Each 
endpiece  also  allowed  for  the  use  of  a splitter  plate  in  various  posi- 


tions. 
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HI . Description  of  Apparatus 


Wind  Tunnel 

Wind  tunnel  tests  were  conducted  in  the  A FIT  Five- Foot  Wind 


Tunnel.  The  tunnel,  an  open  return  induction  type,  has  a maximum 
flow  speed  of  250  miles  per  hour.  Two  large  plyboard  side  panels 
were  installed  in  the  five-foot  diameter  test  section  to  simulate  a two- 
dimensional  test  section  60  in.  high  and  30  in.  wide.  The  trailing 
edges  of  the  side  panels  were  hinged  flaps  and  attached  to  servos  to 
control  the  angle.  Four  pitot-  static  tubes  were  installed  at  the 
entrance  of  the  test  section  on  the  bottom,  top,  right  and  left  sides, 
and  measured  the  local  dynamic  pressure  q.  To  insure  that  these 
four  local  q's  remained  essentially  equal,  the  hinged  flaps  of  the  side 
panels  were  adjusted  with  the  servos.  As  a result,  a more  uniform 

flow  was  maintained  in  the  test  section. 

Eight  static  ports  evenly  spaced  around  the  circumference  of 
the  tunnel  measured  the  static  pressure  at  the  entrance  to  the  test 
section.  The  difference  between  the  static  pressure  and  the  atmos- 
pheric pressure  was  recorded  on  a micro-manometer  filled  with 
water  and  designated  the  tunnel  head  or  "tunnel  q.  n 


Airfoil 


The  airfoil  model  was  a 20  percent  thick,  five  percent 
cambered,  elliptical  airfoil.  The  model  had  a span  of  2.  17  feet  and 
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a chord  of  1.67  feet.  A sketch  of  the  cross  section  of  the  airfoil  is 
shown  in  Fig.  1. 

The  airfoil  was  equipped  with  48  surface  static  pressure  taps 
used  to  measure  the  pressure  distribution  around  the  airfoil.  Forty- 
four  of  the  taps  were  located  along  center  span  and  concentrated 
toward  the  leading  and  trailing  edges  where  greater  pressure  changes 
occur.  Four  taps  were  located  off  the  center  of  the  span  near  the 
blowing  slot  to  check  the  spar.wise  uniformity  of  the  flow.  Three 
pitot  tubes  were  located  inside  the  plenum  chamber  of  the  airfoil  to 
measure  plenum  total  pressure  and  to  check  the  spanwise  uniformity 
of  the  flow  exiting  the  slot. 

Large  circular  aluminum  sideplates  were  fitted  to  each  side  of 

the  airfoil  model.  They  were  — of  an  inch  thick  and  beveled  to  30 

16 

degrees  at  the  edges.  The  purpose  of  these  sideplates  was  to  strip 
off  the  boundary  layer  that  formed  on  the  plyboard  side  panels  and  to 
aid  in  making  the  flow  more  nearly  two-dur  ensional.  The  left  side- 
plate  also  served  as  an  attachment  point  foi  <-L  :es  from  the  angle-of- 
attack  drive  mechanism.  The  cables  were  attached  to  the  front  and 
rear  edges  of  the  sideplate  and  were  run  through  the  tunnel  floor  to  a 
motor  driven  gear  box.  The  gear  box,  in  turn,  was  connected  to  a 
revolution  counter  calibrated  to  an  accuracy  of  three  minutes  of  arc. 

The  blowing  or  secondary  air  system  in  the  model  consisted  of 
a fiberglass  plenum  chamber  and  annealed  copper  pipe.  Both 
chamber  and  pipe  extended  the  span  of  the  model.  A cross  section  of 
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the  plenum  chamber  and  inlet  pipe  also  appear  in  Fig,  1.  Secondary 
air  entered  through  ihc  left  side  of  the  model  into  the  inlet  pipe. 

From  the  pipe,  the  airflow  entered  the  plenum  chamber,  designed  as 
a diverging-converging  nozzle.  The  minimum  area  of  tie  converging 
portion  of  the  chamber  served  as  the  blowing  slot  from  which  air 
exited  onto  the  trailing  edge  surface. 

Four  endpieces,  constructed  to  be  mounted  on  the  airfoil,  were 
tested.  The  endpieces  were  designated  as  Endpieces  A,  B,  C,  and  D. 
The  different  endpieces  appear  in  Figures  2 and  3 and  are  listed  in 
Table  1. 


T a bl  e 1 

The  Four  Endpieces  Used  in  the  Wind  Tunnel  Tests 


Endpiece 


Slot  Position 

Slot  Angle 

Contour 

90%  chord 

+ 5 degrees 

elliptical 

90%  chord 

- 33  degrees 

ell  iptical 

96%  chord 

t 5 degrees 

c i rcula  r 

97%  chord 

+ 5 degrees 

ell  iptical 

Endpiece  A was  similar  in  construction  to  the  trailing  edge  configura- 
tion used  by  Stevenson  and  Rhynard.  The  slot  angle  on  Endpiece  13 
was  constructed  to  be  more  nearly  tangential  to  the  surface.  The 
circular  contour  on  Endpiece  C had  a radius  of  0.  95  inches.  The 
slots  in  the  four  endpieces  were  located  on  the  upper  surface  and 


extended  across  the  span  of  the  endpieces.  Each  slot  was  0.  020 
inches  high,  fixed  at  a height  to  chord  ratio  of  0.001,  and  uniform  to 
+_  0.0015  inches  along  the  span. 

Three  configurations  of  splitter  plates  were  tested  on  each  end- 
piece.  They  were  designated  as  No  Splitter  Plate,  Aft  Splitter  Plate, 
and  Forward  Splitter  Plate.  The  splitter  plate  used  had  a chord 
length  of  1.  5 inches.  The  No  Splitter  Plate  configuration  had  no 
splitter  plate  on  the  lower  surface  and  meant  that  the  airfoil  was 
clean.  The  Aft  Splitter  Plate  configuration  had  the  splitter  plate  at 
the  99  percent  chord  position  and  at  an  angle  of  -45  degrees  from  the 
chord  line.  The  splitter  plate  in  the  I-orward  Splitter  Plate  config- 
uration was  located  at  the  95.  3 percent  chord  position  and  at  angle  of 
-60  degrees  from  the  chord  line.  The  latter  two  configurations  were 
chosen  because  they  were  found  by  Rhynard  (Ref  4)  to  be  optimum 
configurations.  These  configurations  appear  in  Figures  2 and  3. 


Flowmeter 

The  flowmeter  used  to  measure  the  mass  flow  rate  was  a one 


1 


inch  diameter  tube  venturi  with  a throat  diameter  of  0.  50  inches.  It 
was  calibrated  against  a National  Bureau  of  Standards  venturi.  The 
temperature  of  the  air  flowing  through  the  venturi  was  measured  by 
a copper-constantan  thermocouple  inserted  just  upstream  of  the 
throat. 
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W ake  Survey  Rake 

A total  head  wake  survey  rake,  designed  to  measure  the 
momentum  deficit  of  the  airfoil  wake,  was  positioned  37  inches,  or 
1.85  chord  lengths,  behind  the  airfoil.  The  rake  had  96  total  head 
tubes  and  two  static  tubes.  The  tubes  were  0.  0625  inches  in  outside 
diameter  and  spaced  0.  25  inches  apart.  The  airfoil  section  of  the 
rake  spanned  the  tunnel  from  top  to  bottom  and  was  adjustable  within 
certain  limitations. 


Manometers 


Two  manometer  banks  were  used  to  measure  the  pressure 
distributions  on  the  airfoil  surface  and  the  wake  survey  rake.  A 96 
inch,  100-tube  bank  of  vertical  manometers  containing  alcohol  was 
connected  to  the  pressure  taps  on  the  airfoil.  A 30  inch,  100-tube 
manometer  bank  containing  red  oil  was  connected  to  the  wake  survey 
rake.  The  bank  was  inclined  60  degrees  from  the  vertical  so  that 
pressure  changes  on  the  rake  could  be  read  more  accurately. 

Two  60  inch  vertical  mercury  manometers  were  used  with  the 
flowmeter  to  measure  the  secondary  air  flow  rate.  Finally,  a 
mercury  manometer  bank  was  connected  to  the  plenum  chamber  pitot 
tubes  to  measure  the  total  pressure  in  the  chamber. 
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Pitot  Tube 

A small  pitot  tube  was  constructed  to  measure  the  slot  span- 
total  pressure  distribution.  The  outside  diameter  was  0.  018 


wise 


inches  s 


o that  the  pitot  tube  was  able  to  fit  into  the  slot. 
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IV.  Experimental  Procedure 

The  same  general  test  procedure  was  ued  on  each  of  the  four 
model  and  endpiece  configurations.  Once  the  endpiece  was  put  on  the 
model  it  was  checked  for  leaks  and  uniformity  of  spanwise  flow.  The 
tunnel  speed  was  then  brought  up  to  1.  25  inches  of  H^O,  correspond- 
ing to  a nominal  74  feet  per  second.  The  geometric  angle  of  attack 
was  varied  from  -6  to  +6  degrees  in  increments  of  two  degrees.  At 
each  angle  of  attack  the  two  manometer  banks  were  photographed  and 
data  recorded  from  the  flowmeter  and  plenum  chamber  manometers. 
After  the  angle  of  attack  sequence  was  completed,  the  secondary  air 
system  was  started  and  blowing  rate  adjusted  to  22  inches  of  mercury 
on  the  flowmeter  manometer.  This  corresponded  to  a blowing 
momentum  coefficient,  CM,  of  approximately  0.03.  The  angle  of 
attack  sequence  was  then  repeated  for  this  blowing  rate.  The  tunnel 
speed  was  reduced  to  zero  and  the  splitter  plate  was  mounted  in  the 
aft  position.  The  complete  process  was  repeated  for  this  configura- 
tion and  then,  after  the  splitter  plate  was  remounted  in  the  forward 
position,  was  repeated  again  for  that  configuration.  Test  runs  were 
repeated  for  Endpieces  A and  B as  a check  for  accuracy. 
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V,  Data  Reduction 

Data  obtained  in  the  wind  tunnel  tests  were  reduced  to  three 
primary  parameters.  The  parameters  were  the  blowing  momentum 
coefficient,  the  section  lift  coefficient,  and  the  section  total  drag 
coefficient.  These  parameters  were  then  used  to  calculate  the  liit- 
to-drag  ratios  of  the  various  model  and  endpiece  configurations. 

Momentum  Coefficient 

rlhe  momentum  coefficient  is  a measure  of  the  amount  of  blow- 
ing applied  to  a circulation  controlled  airfoil.  It  is  defined  as 


V 


q0c 


(i) 


where  m is  the  mass  flow  rate  of  air  per  unit  span,  Vj  is  the  velocity 
of  the  air  jet  at  the  slot,  qQ  is  the  free  stream  dynamic  pressure,  and 
c is  the  chord  length  of  the  airfoil  (Ref  5 : 1 6 ) . 


Section  Lift  Coefficient 


The  section  lift  coefficient  was  calculated  by  the  equation: 


Cj  = Cn  cos  a 


(2) 


where  a is  the  geometric  angle  of  attack  and  Cn  is  the  section  normal 
force  coefficient.  Cn  was  calculated  by  numerical  integration  of  the 
pressure  coefficients  around  the  airfoil.  The  pressure  coefficient  is 

14 


GAE/AE/75D-12 


defined  as 


C = 


?!  - Po 


Integration  was  accomplished  according  tc  the  equation: 


c-  /'(c"  ■ c-)  ‘ (?) 


where  C and  C_  are  pressure  coefficients  on  the  lower  and  upper 
Pi  Pu 


surfaces,  respectively,  and  — is  the  distance  on  the  chord  line 

c 


(Ref  5:17).  This  was  performed  on  the  Hewlett-Packard  9100A 


Calculator  with  9107A  Digitizer. 


Section  Total  Drag  Coefficient 


The  section  total  drag  coefficient  takes  into  account  the  section 


profile  drag  coefficient,  an  intake  penalty  coefficient,  and  the  blowing 


momentum  coefficient.  The  section  total  drag  coefficient  is  deter- 


mined by  the  equation: 


Cdfc  - cd  + + 


m V0 


The  section  profile  drag  coefficient  was  calculated  from  the  momen- 


tum method  and  compares  the  momentum  ahead  of  the  model  with  the 


momentum  behind  the  model.  It  is  defined  as 


Cd  = t/ 


h/  p:__CL 

VPo 
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where  q and  qQ  are  the  wake  and  free  stream  dynamic  pressures, 
respectively,  and  dy  is  the  incremental  distance  between  tubes  on  the 
rake  (Ref  5:18).  The  integration  was  performed  by  the  trapezoidal 
rule  of  integration  on  the  Hewlett-Packard  9100A  Calculator  with 
9107A  Digitizer.  An  intake  penalty  coefficient,  from  Englar  (Ref 
1:10),  is  applied  because  there  is  an  increase  in  drag  associated  with 
the  production  of  the  blowing  air  from  the  free  stream  through  the 
system  to  the  slot.  The  momentum  coefficient  term,  , is  added 
because  some  of  the  momentum  originated  inside  the  model  and  this 


momentum  actually  decreased  the  momentum  deficit  that  represents 
the  drag.  The  two  additional  terms  were  included  in  calculation  of 
the  section  total  drag  coefficient  in  order  to  allow  for  comparison  to 
unblown  airfoils. 
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stream  atmosphere  accurately.  The  presence  of  the  model  in  the 
test  section  caused  the  free  stream  velocity  to  increase  as  it  flowed 
over  the  model  and  the  presence  of  the  floor  and  ceiling  prevented 
normal  curvature  of  the  flow.  To  account  for  these  characteristics, 
streamline  curvature  corrections  were  applied  to  C\  while  wake  and 
solid  blocking  corrections  were  applied  to  Cp  Cd,  qQ,  and  Re.  A 
wake  survey  rake  correction  factor  w?s  checked  to  apply  to  the  pres 
sure  readings  but  it  was  negligible. 
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VI.  Analytic  Stud-\ 


The  analytic  study  consisted  of  the  development  and  manipula- 
tior  of  a potential  flow  computer  program.  Its  purpose  was  to  pre- 
dict results  obtained  in  wind  tunnel  testing  of  the  four  endpieces,  with 
various  splitter  plate  configurations.  It  was  based  on  potential 
theory  and  assumed  inviscid,  incompressible  flow. 


Development 

The  airfoil  surface  was  segmented  into  60  panels  and  panel  size 
became  smaller  at  the  leading  and  trailing  edges  so  as  to  weight  these 
more  sensitive  areas.  A vortex  was  located  on  each  panel  at  the  ^ 

3 

chord  point  and  a control  point  was  located  at  the  — chord  point.  This 
paneling  scheme  modeled  the  airfoil  configuration  without  a splitter 
plate.  To  model  the  airfoil  with  a splitter  plate,  two  additional  panels 
were  affixed  to  the  lower  surface  in  the  appropriate  position  near  the 
trailing  edge.  The  three  configurations  of  no  splitter  plate,  aft  split- 
ter plate,  and  forward  splitter  plate  are  represented  in  Fig.  4.  The 
splitter  plates  in  the  program  are  not  the  same  size  as  in  the  tests 
because  it  was  found  that  the  thickness  and  length  do  not  have  a 
significant  effect  in  potential  flow  theory. 

The  next  step  in  the  development  was  to  solve  for  the  particular 
vortex  strengths,  with  the  boundary  condition  that  normal  velocity 
through  the  surface  was  zero.  From  these,  the  lift  coefficients 
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c 


C 


were  c 


alculated.  The  potential  function  for  a vortex  is  given  by  the 


equation: 


<t>  = 

^ v 


tan 


1 r Xg  ~ Xvi 

[ye  - YvJ 


(8) 


where  x and  yc  are  the  coordinates  of  the  control  point  and  xv  and  yv 


are 


the  coordinates  of  the  vortex  influencing  that  control  point.  From 
the  potential  function  the  u and  v components  of  velocity  relative  to 
the  chord  line  are  given  by 

(yc  - Yv) 


r 

U 2* 


(yc  - Yv)2  + (XC  - xv)' 


(9) 


v = 


-r 


2* 


l 


(xc  - Xv) 


(yc  - yv)  + (xc  - 


(10) 


The  tangential  and  normal  velocities  on  the  surface  of  the  air- 
foil due  to  the  vortices  are  then  defined  as 


VT  = u cos  0 + v sin  0 


VN  = u sin  0 - v cos  0 


(ID 

(12) 


where  0 is  the  angle  of  the  surface,  relative  to  the  chord  line,  at  the 
control  point.  The  normal  velocity  is  also  equal  to  the  negative  of 
the  free  stream  velocity  component.  To  account  for  various  angles 
of  attack  this  equation  becomes 

VN  = - UQ  sin  ( 0 - a ) (13) 


i 
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where  Uq  is  the  free  stream  velocity  and  a is  the  geometric  angle  of 
attack. 

Sixty  panels  yielded  60  linear  equations  such  that 

[CyN]  [r]  = - u0  [sin  (6-  a)]  (14) 

where  Cy^  is  the  coefficient  of  normal  velocity  matrix,  r is  the 
vortex  strength  matrix,  and  sin  ( 9 - a ) is  the  free  stream  component 
of  velocity  matrix.  Solving  this  system  of  linear  equations  simul- 
taneously yielded  the  vortex  strengths.  The  vortex  strengths,  T,  were 
then  substituted  into  Eqs  (9)  and  (10)  and  the  results  were  substituted 
into  Eq  (11).  These  equations  gave  the  tangential  velocity  at  each 
control  point  induced  by  the  vortices  at  each  panel.  The  V^'s  were 
then  related  to  the  dimensionless  pressure  coefficients  by  the 
equation: 


In  order  to  obtain  dimensionless  lift  coefficients  the  Cp's  were  then 
integrated  around  the  airfoil  for  each  angle  of  attack. 

Prediction  Technique 

The  prediction  technique  developed  was  limited  to  this  study. 

It  was  based  on  initial  wind  tunnel  results  of  one  airfoil  and  endpiece 
configuration  and  attempted  to  predict  the  results  of  using  the  other 
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three  endpieces.  This  was  accomplished  by  matching  the  lift  coef- 
ficient results  of  Endpiece  A and  making  changes  in  the  program  to 


simulate  Endpieces  B,  C,  and  D.  To  match  the  lift  coefficient  test 
results  of  Endpiece  A in  the  program,  the  on  the  upper  surface 
was  controlled  and  thus  the  circulation,  also.  First,  a particular 
panel  control  point  was  chosen  representing  the  location  of  the  actual 
blowing  slot.  Corresponding  to  this  point  the  one  particular  linear 
equation  of  normal  velocity  in  the  matrix  was  rewritten  as  tangential 
velocity.  It  was  set  equal  to  a variable  by  the  equation: 


Yry  = A 


where  A appears  in  the  sin  (0  - a)  matrix.  Then  the  vortex  for  that 

particular  panel  was  relocated  below  the  surface,  perpendicular  to  the 

control  point,  and  — of  a panel  length  away  from  it.  This  was  done  so 

4 

the  vortex  would  have  more  influence  on  that  control  point.  Lift  coef- 
ficients at  angles  of  attack  for  Endpiece  A were  then  obtained  by  using 
certain  values  of  A.  For  a fixed  angle  of  attack,  the  more  negative  A 
becomes,  the  greater  the  V-j-.  Inconsistencies  in  lift  curve  slopes 
between  test  data  and  the  programs  necessitated  changing  the  values  of 
A throughout  a range  of  angles  of  attack,  for  most  cases.  These 
values  of  A appear  in  Table  II.  Thus,  two  sets  of  A values  were 


obtained  corresponding  to  the  lift  curve  slopes  of  Endpiece  A with  no 
blowing  and  blowing  at  a of  0.  03, 
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The  next  step  was  to  predict  the  test  results  with  Endpieces  B, 

C,  and  D.  To  accomplish  this,  the  program  haa  to  incorporate  the 
changes  of  a new  slot  angle,  a different  trailing  edge  contour,  and  a 

{] 

second  slot  position. 

The  program  did  not  have  an.  inherent  capability  of  a slot  angle 
change.  To  simulate  the  more  tangential  blowing  of  Endpiece  B,  the 
values  of  A were  made  more  negative  and  thus,  V<p  was  increased. 

This  was  based  on  the  assumption  that  the  more  tangential  the  blowing, 
the  greater  the  velocity  immediately  downstream  of  the  slot.  These 
additional  values  of  A also  appear  in  Table  II. 

The  circular  contour  of  Endpiece  C was  put  into  the  program  by 

!i 

changing  the  paneling  scheme  on  the  trailing  edge.  Instead  of  an  equa- 
tion of  an  ellipse,  an  equation  of  a circle  was  used  to  define  the 
surface  for  the  last  few  panels.  The  A values  obtained  for  Endpiece  A 
were  then  used  in  this  case  to  specify  V<p  at  the  slot  location.  These 
were  used  assuming  that  the  amount  of  blowing  was  essentially  equal 
for  both  cases . 

The  slot  position  change  to  97  percent  of  chord  in  Endpiece  D 
was  accomplished  by  choosing  a different  panel  and  repeating  the 
previous  procedure  for  Endpiece  A.  A control  point  and  vortex  loca- 
tion were  chosen  one  panel  downstream  of  the  one  used  for  Endpiece 
A.  Here  again,  to  achieve  lift  curve  slopes,  the  same  A values 
obtained  for  Endpiece  A were  used  to  simulate  the  comparable 
amounts  of  blowing. 
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VII.  Results  and  Discussion 


General  Observations 

Several  phenomena  were  observed  in  the  wind  tunnel  testing  of 
the  model  and  endpiece  configurations.  The  lift  curves  of  the  four 
endpieces  without  blowing  were  essentially  equ  1.  This  similarity  in 
lift  curves  appears  in  Figures  5 through  8 and  shows  that  the  config- 
uration changes  of  each  endpiece  did  not  have  a significant  effect  on  the 
flow  without  blowing.  However,  with  blowing,  Figures  9 through  12, 
the  endpieces  had  different  lift  curves  showing  the  effect  of  the  config- 
uration changes.  Also,  in  these  curves  a tendency  of  flow'  separation 
at  higher  angles  of  attack  was  observed.  This  stall  regime  generally 
occurred  from  +2  degrees  to  +5  degrees  geometric  angle  of  attack. 

The  drag  curves,  Figures  13  through  20,  display  a tendency  of 
increasing  drag  towards  extreme  negative  and  positive  angles  of 
attack.  These  curves  show  that  the  smallest  amounts  of  drag  occurred 
around  zero  angle  of  attack  and  the  largest  amounts  occurred  towards 
+6  degrees  angle  of  attack. 

The  lift-to-drag  ratio  curves.  Figures  21  through  28,  exhibit  a 
tendency  of  decreasing  L/ D at  the  extreme  negative  and  positive  angles 
of  attack.  This  was  primarily  due  to  the  increasing  amounts  of  drag 
in  these  areas. 
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Lift  Results 


The  lift  results  are  presented  in  Figures  5 through  12  as  section 
lift  coefficient  plotted  versus  geometric  angle  of  attack.  The  first  four 
are  for  the  no  blowing  case  on  each  endpiece  and  the  last  four  are  for 
the  blowing  case  on  the  same  four  endpieces. 

In  the  no  blowing  cases,  the  Cj's  of  the  forward  splitter  plate 
configuration  increased,  on  the  average,  by  1.  14  over  the  no  splitter 
plate  configuration.  The  aft  splitter  plate  configuration's  Cj's 
increased,  on  the  average,  by  0 87  over  the  no  splitter  plate  config- 
uration. Thus,  a higher  section  lift  coefficient  was  achieved  by  use  of 
the  forward  splitter  plate  on  each  endpiece  without  blowing. 

In  the  blowing  cases  the  results  were  much  the  same.  For  the 
forward  splitter  plate,  the  Cj's  increased,  on  the  average,  by  0.78 
from  the  no  splitter  plate  values.  And  the  aft  splitter  plate  Cj's 
increased,  on  the  average,  by  only  0.  69  from  the  no  splitter  plate 
values.  Thus,  the  use  of  the  splitter  plate  in  the  forward  position  pro- 
vided the  highest  section  lift  coefficients  throughout  the  tested  range  of 
angles  of  attack.  In  the  no  blowing  case,  these  values  represented  a 
280  percent  average  increase  in  Cj  over  the  no  splitter  plate  configura- 
tion and  a 60  percent  average  increase  in  Cj  over  the  aft  splitter  plate 
configuration.  In  the  blowing  case,  the  increases  over  the  no  splitter 
plate  and  aft  splitter  plate  configurations  were,  on  the  average,  190 
percent  and  20  percent,  respectively.  These  percentages  were 


calculated  for  the  middle  range  angles  of  attack. 
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Drag  Results 

The  drag  results  are  presented  in  Figures  13  through  20  as 
section  total  drag  coefficient  plotted  versus  geometric  angle  of  attack. 
The  first  four  are  without  blowing  and  the  last  four  are  with  blowing. 

In  the  no  blowing  cases,  the  configuration  without  the  splitter 
plate  provided  the  lowest  amounts  of  drag.  The  aft  splitter  plate 
caused  higher  amounts  of  drag  and  the  Cdj.'s  increased  an  average  95 
percent  over  the  no  splitter  plate  case.  The  forward  splitter  plate  had 
even  higher  Cd^'s  and  these  were  an  average  130  percent  higher  than 
the  no  splitter  plate  case.  These  percentages  were  calculated  for  the 
middle  range  angles  of  attack. 

The  drag  results  in  the  blowing  cases  were  somewhat  less  defin- 
itive. In  some  cases  the  aft  splitter  plate  had  the  lowest  Cd^-'s.  This 
occurred  around  the  zero  angle  of  attack.  Figures  17  and  19,  and 
around  the  extreme  negative  angles  of  attack.  Figure  20.  The  no 
splitter  plate  configuration  had  the  least  amount  of  drag  when  blowing 
was  not  used  but  when  blowing  was  applied,  the  losses,  due  to  the 
mixing  of  air  flows  below  the  trailing  edge,  caused  the  drag  to  increase. 
This  increase  was  even  more  than  that  of  the  aft  splitter  plate  config- 
uration because  the  plate  actually  prevented  mixing  of  flows  and 
reduced  these  losses.  However,  Figures  18,  19,  and  20  show  that  at 
most  of  the  argles  of  attack  the  Cd^'s  were  lowest  for  the  no  splitter 
plate  case.  In  all  cases  the  drag  was  highest  for  the  forward  splitter 
plate  configuration. 
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Thus,  the  model  and  endpiece  configurations  achieved  the  lowest 


amounts  of  drag  for  the  most  part  when  no  splitter  plate  was  used  and 


had  the  highest  amounts  when  a splitter  plate  was  used  in  the  forward 


position. 


Lift-to-Drag  Ratio  Results 


The  lift-to-drag  ratio  results  are  presented  in  Figures  21 


through  28  as  L/D  plotted  versus  geometric  angle  of  attack.  The  first 


four  are  without  blowing  and  the  last  four  are  with  blowing. 


In  the  no  blowing  cases,  each  endpiece  had  a slightly  different 


result.  For  Endpiece  A,  Fig.  21,  the  forward  splitter  plate  produced 


the  highest  JL/D's  throughout  the  range  of  angles  of  attack.  This  was 


primarily  due  to  the  fact  that  it  had  much  greater  C^'s  than  the  other 


lower  surface  configurations.  On  Endpiece  B,  Fig.  22,  the  forward 


splitter  plate  produced  high  L/D’s  at  the  negative  angles  of  attack  but 


then  dropped  below  the  higher  L/D  levels  of  the  aft  splitter  plate 


and  no  splitter  plate  configurations.  Up  to  the  stall  region,  the  high- 


est L/D  values  of  Endpiece  C,  Fig.  23,  were  produced  with  the  aft 


splitter  plate.  On  Endpiece  D,  Fig.  24,  the  aft  splitter  plate  had 


slightly  greater  L/D's  at  the  positive  angles  of  attack  but,  due  to  the 


relatively  small  difference  between  these  values  and  the  forward  split- 


ter plate  values,  a conclusive  statement  about  the  relative  merit  of 


either  configuration  cannot  be  made.  Throughout  the  no  blowing  cases. 


the  no  splitter  plate  configuration  had  the  lowest  L/D's,  due  primarily 
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to  low  Cj's.  The  L/D's  were  still  on  the  increase  at  the  extreme  posi- 


tive angles  of  attack  thus  suggesting  more  study  at  even  higher  angles 


of  attack. 


The  lack  of  uniformity  and  definite  trends  in  the  data  for  the  no 


blowing  case  can  be  attributed  to  significant  sensitivity  in  the  drag  data 


reduction  technique. 


In  the  blowing  cases,  however,  the  results  are  more  conclusive. 


On  Endpieces  A,  C,  and  D,  Figures  25,  27,  and  28,  the  use  of  the  aft 


splitter  plate  produced  the  highest  L/ D values.  Such  results  were  due 


to  the  fact  that  the  plate  in  the  aft  position  significantly  reduced  the 


mixing  losses  enough  to  reduce  the  values.  On  Endpiece  B,  Figur 
26,  the  forward  splitter  plate  actually  had  slightly  greater  L/D's  than 


the  aft  splitter  plate  but  the  relative  difference  was  not  enough  to  make 


a definite  statement  about  this  configuration. 


Thus,  throughout  the  tested  range  of  angles  of  attack,  the  highest 


lift-to-drag  ratios  were  achieved  on  the  model  and  endpiece  configura- 


tions with  the  splitter  plate  in  the  aft  position.  In  the  middle  range 


angles  of  attack  these  L/D's  of  the  aft  splitter  plate  configuration 


represented  an  average  increase  of  37  percent  over  the  forward  splitter 


plate  configuration  and  an  average  increase  of  64  percent  over  the  no 


splitter  plate  configuration.  This  was  true  in  the  blowing  cases  and  in 


some  of  the  no  blowing  cases.  In  only  one  case  did  the  forward  split- 


ter plate  or  no  splitter  plate  configurations  exhibit  significantly  higher 


L/D  values  than  the  aft  splitter  plate  values. 
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Thus,  whether  the  airfoil  had  the  slot  located  at  the  96  or  97 
percent  chord  position,  at  an  angle  of  +5  or  -33  degrees,  or  had  a 
circular  or  elliptical  contour,  the  splitter  plate  in  the  aft  position 
provided  the  maximum  lift-to-drag  ratios.  The  effect  of  blowing  was 
to  generally  increase  the  L/D  values  over  those  of  the  no  blowing 
cases.  This  was  true  with  the  exception  of  a few  cases  where  the 
secondary  flow  was  non-uniform  through  the  slot.  Further  discus- 
sion of  the  effects  of  blowing  can  be  found  in  the  investigation  by 
Oxford  (Ref  3). 


Potential  Flow  Results 

The  results  of  the  potential  flow  computer  program  are  pre- 
sented in  Figures  29  through  34.  The  first  three  show  the  results  of 
the  lift  curve  matching  routine  and  the  last  three  show  the  results  of 
the  prediction  technique. 

The  matching  routine,  accomplished  with  the  values  of  A in 
Table  II,  assumed  that  potential  curves  would  not  exhibit  the  decrease 
in  Cj  that  normally  occurs  when  the  flow  separates.  The  values  used 
were  obtained  with  a linear  equation  and  the  results  show  somewhat 
linear  potential  flow  curves  that  follow  the  test  lift  curves  closely,  up 
to  the  stall  regions. 

The  prediction  technique  for  each  splitter  plate  configuration 
was  based  on  the  A values  established  for  Endpiece  A,  with  blowing. 
Figures  32,  33,  and  34  show  that  when  the  contour  was  varied  from 
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Erdpiece  A to  C,  the  test  lift  curves  changed  appreciably  but  the 
potential  flow  curves  did  not.  This  result  was  possibly  due  to  the  fact 
that  the  flow  in  potential  theory  is  defined  to  be  inviscid.  Under  real 
flow  conditions  a circular  contour  makes  the  blowing  and  subsequent 
reattachment  more  effective  and  thus  increases  the  Cq's,  but  in 
potential  flow  a contour  change  has  an  insignificant  effect. 

When  the  slot  position  was  moved  from  96  to  97  percent  of  the 
chord,  Endpieces  A to  D,  the  real  flow  test  Cq's  increased  appreciably. 
The  potential  flow  Cq's  did  not  increase  but  actually  decreased  at  nega- 
tive angles  of  attack.  This  was  due  to  the  fact  that  the  potential  flow 
program  did  not  accurately  simulate  real  flow  conditions  at  and  around 
the  splitter  plates  and  lower  surface  stagnation  points. 

The  results  of  the  prediction  technique  show  that  it  is  not  useful 
for  prediction  purposes  because  it  is  constrained  to  this  study  and  is 
based  on  initial  wind  tunnel  data.  A true  prediction  technique  would 
not  be  based  on  initial  data  and  would  be  capable  of  broader  usage. 
Thus,  the  results  imply  that  the  program  is  more  justifiable  as  a 
matching  routine  as  it  can  match  any  of  the  lift  coefficient  test  results. 
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VIII.  Conclusions 


A two-dimensional  wind  tunnel  and  analytic  study  to  determine 


the  effects  of  splitter  plate  position  on  trailing  edge  modifications  of 


a circulation  controlled  airfoil  resulted  in  the  following  conclusions. 


1.  The  section  lift  coefficient  increases  as  the  splitter  plate  is 


moved  forward  toward  the  95.  3 percent  chord  position. 


2.  The  section  total  drag  coefficient  increases  as  the  splitter 


plate  is  moved  forward  toward  the  95.  3 percent  chord 


position. 


3.  The  lift-to-drag  ratio  is  a maximum  when  the  splitter  plate 


is  located  at  the  99  percent  chord  position. 


4.  The  above  conclusions  are  valid  for  a circulation  controlled 


airfoil  with  variations  in  slot  position,  slot  angle,  and  trail- 


ing edge  contour. 


5.  Computer  methods  based  on  potential  flow  theory  can  be 


used  to  match  lift  coefficients  obtained  in  wind  tunnel  test- 


ing of  a circulation  controlled  airfoil. 


Computer  methods  based  solely  on  potential  flow  theory 


cannot  be  used  to  predict  lift  coefficient  results  obtained  in 


wind  tunnel  testing  of  a circulation  controlled  airfoil. 
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IX.  Recommendations 


It  is  recommended  that  further  wind  tunnel  tests  and  analytic 
studies  be  made  on  circulation  controlled  airfoils.  Wind  tunnel  test- 


ing should  include: 

1.  A study  of  the  effects  of  a splitter  plate  that  is  free  to 
rotate  and  seek  its  own  optimal  angle  with  the  chord  line. 

2.  A study  of  the  effects  of  such  a splitter  plate  configuration 
on  further  modifications  of  the  trailing  edge  of  the  airfoil. 

Further  analytic  investigation  should  include  a potential  flow  study  that 
incorporates  the  effects  of  viscosity. 
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Fig.  6 The  Change  in  vs  a Due  to  the  Different  Splitter  Plate 
Configurations  on  Endpiece  B,  Without  Blowing 
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Fig.  12  The  Change  in  C1  vs  a Due  to  the  Different  Splitter  Plate 
Configurations  on  Endpiece  D,  With  Blowing 
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Fig.  14  The  Change  in  vs  a Due  to  the  Different  Splitter  Plate 
Configurations  on  Endpiece  B,  Without  Blowing 
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Fig.  15  The  Change  in  Cd(.  vs  a Due  to  the  Different  Splitter  Plate 
Configurations  on  Knclpiece  C,  Without  Blowing 
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Fig.  22  The  Change  in  L/D  vs  a Due  to  the  Different  Splitter  Plate 
Configurations  on  Endpiece  B,  Without  Blowing 


LI FT-TO-DRRG  RATIO  { L/D) 


GAE/AE/75D-12 


ENDPIECE  D 

O-NO  SPLITTER  PLATE 
CD-flFT  SPLITTER  PLATE 
O-FHD  SPLITTER  PLATE 
Re  = 7.36  X 10  6 
Cu  = 0.0 


60. 


40.H 


-10.  J 


ANGLE  OF  ATTACK  (DEG) 


Fig.  24  The  Change  in  L/D  vs  a Due  to  the  Dili  rent  Splitter  Plate 
Configurations  on  Endpiece  D,  Without  Blowing 
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A Comparison  of  the  Results  of  Configuration  Changes  Made 
in  the  Potential  Flow  Program  and  in  the  Wind  Tunnel  Tests, 
With  Aft  Splitter  Plate 
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Table  II 

The  Values  of  the  Potential  Flow  Program 
Variable  A Used  to  Match  Lift  Curve  Slopes 


Angles  of  Attack 


mm 

-4° 

-2° 

0° 

wm 

Esa 

6° 

Endpiece  A,  = 0.  0 

No  Splitter  Plate 

0.  85 

EH 

0.  35 

0.  10 

-0.  15 

-0.  40 

-0.  65 

Aft  Splitter  Plate 

-0.  82 

-1.18 

-1.  36 

-1.  54 

-1.  72 

-1.  90 

Forward  Splitter  Plate 

-1.12 

m 

-1.  52 

-1.  72 

-1.  92 

-2.  12 

-2.  32 

Endpiece  A,  C ^ = 0.  03 

No  Splitter  Plate 

-0.  85 

-0.  85 

-0.  85 

-0.  85 

-0.  85 

-0.  85 

Aft  Splitter  Plate 

i 

eo 

O 

-2.  30 

-2.  20 

-2.  10 

-2.  00 

-1.  90 

-1 . 80 

Forward  Splitter  Plate 

-2.  68 

-2.  53 

-2.  38 

-2.  23 

-2.  08 

-1.93 

-1 . 78 

Endpiece  B,  = 0.  03 

No  Splitter  Plate 

-0.  70 

-0.  85 

-1.  00 

-1.15 

-1.  30 

-1. 45 

-1.  60 

Aft  Splitter  Plate 

-2.  40 

-2.  40 

-2.  40 

-2.  40 

-2.  40 

-2.  40 

-2.  40 

Forward  Splitter  Plate 

-2.  75 

-2.  75 

-2.  75 

-2.  75 

-2.  75 

-2.  75 

-2.  75 
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